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ABSTRACT 



Context. Many binary stellar systems in which the primary star is beyond the asymptotic giant branch (AGB) evolutionary phase 
show significant orbital eccentricities whereas current binary interaction models predict their orbits to be circularised. 
Aims. In the search for a mechanism to counteract the circularising effect of tidal interaction we analyse how the orbital parameters 
in a system are modified under mass loss and mass exchange among its binary components. 

Methods. We propose a model for enhanced mass-loss from the AGB star due to tidal interaction with its companion, which allows a 
smooth transition between the wind and Roche-lobe overflow mass-loss regimes. We explicitly follow its effect along the orbit on the 
change of eccentricity and orbital semi-major axis, as well as the effect of accretion by the companion. We calculate timescales for 
the variation of these orbital parameters and compare them to the tidal circularisation timescale. 

Results. We find that in many cases, due to the enhanced mass loss of the AGB component at orbital phases closer to the periastron, the 
net eccentricity growth rate in one orbit is comparable to the rate of tidal circularisation. We show that with this eccentricity enhancing 
mechanism it is possible to reproduce the orbital period and eccentricity of Sirius system, which under the standard assumptions of 
binary interaction is expected to be circularised. We also show that this mechanism may provide an explanation for the eccentricities 
of most barium star systems, which are expected to be circularised due to tidal dissipation. 

Conclusions. By proposing a tidally enhanced model of mass loss from AGB stars we find a mechanism which efficiently works 
against the tidal circularisation of the orbit. This mechanism can explain the significant eccentricities observed in binary systems 
containing a white dwarf and a less evolved companion, which are predicted to be circularised due to their proximity, such as Sirius 
and systems with barium stars. 
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1. Introduction 

Detached binary systems containing a white dwarf and a rel- 
atively unevolved companion, i.e., a main-sequence star or a 
(sub)giant, are a useful tool to understand the binary evolution 
of systems with an asymptotic giant branch (AGB) star, given 
that their orbital and chemical properties do not change signifi- 
cantly from the moment that the primary finished its AGB evo- 
lution and became the current white dwarf. An example of such 
a system is Sirius, a 2.1 main sequence star wit h a 1 .05 Mp 
white dwarf companion in a 50-ye ar orbit (e.g., ivan den BosI 
ll960l:lGatewood & Gatewoodlll978[) . Under the standard picture 
of binary evolution with tidal interaction this system should have 
circularised when the primary became an AGB star, as we show 
in ^4.1l of this paper. However, this system has an eccentricity 
e = 0.59 and is not an exception. A similar problem is faced 
when considering the barium-star systems, which are red giants 
with over-abundances of i-process elements (prominently bar- 
ium) with white dwarf companions. The best current explanation 
for the enhancement of .^-process elements in these stars is that 
they accreted mass from their companion when it was an AGB 
star. Due to the large size of AGB stars these systems are ex- 
pected to be circularised by tidal interaction for periods smaller 
than about 3500-4000 days ( Pols et al. 2003.) . However, barium- 
star systems with period as short as 600 days are observed to 
be significantly eccentric ( Jorissen et al.lll998l) . The problems 
stated above are indications that a mechanism must exist that 
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counteracts the circularising effect of the tides by the time the 
primary star is on the AGB. IVan Winckel et al.] ( Il995h propose 
that enhanced mass loss at periastron could enhance the eccen- 
tricity and Soker (2000) shows that this works when the AGB 
star fills its Roche lobe during periastron passages. It has also 
been proposed that the tidal interaction between a binary system 
and a circumbinary disk could account fo r an eccentricity en - 
hancing mechanism ( Waelkens et al.|[r996l : IWaters et al.lll998l) . 
However, Soker (2000) argues that the masses of observed cir- 
cumbinary disks are too small to counteract the circularisation. 

Mass loss from stars in binary systems is commonly treated 
as single-star wind mass-loss as long as the stars are inside their 
Roche lobes, while when one star fills its Roche lobe the mass- 
loss rate increases abruptly to the high values that correspond to 
Roche-lobe overflow. This approximation is fairly accurate for 
stars with a steep density gradient in their atmospheres, but it is 
not appropriate for the case of AGB stars, given their large atmo- 
spheric pressure scale height and their weakly bound envelope. 
The fact that AGB stars undergo dynamical pulsations further re- 
duces density gradient in the layers above the photosphere (e.g., 
Bowen 1988). We propose a prescription for enhanced mass loss 
which smoothly grows from the single star wind mass loss rate 
to the Roche lobe overflow mass loss rate as the radius of the 
star approaches its Roche lobe radius. This gives a variable mass 
loss rate along eccentric orbits which is higher at orbital phases 
closer to the periastron, even if the star does not fill its Roche 
lobe. This effect works in the same way as discussed above, but 
no filling of the Roche lobe is needed, so that (during the AGB 
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phase) this effect is permanently competing against tidal circu- 

lari sation. 

lBoffin& Jorissenl(ll988 l) carried out calculations of the vari- 
ation of orbital parameters considering instantaneous mass trans- 
fer and only linear momentum conservation. Angular momen- 
tum conservation was not taken into account in this exploratory 
study. An extension including angular momentum conservation 
was made by Liu et al. (2000), but they make use of the orbital 
average of the distance between the components and orbital an- 
gular velocity to carry out their calculations, which in the case 
of a variable mass-loss along the orbit does not give the correct 
results. 

In ^ we revise the variation of orbital parameters due to 
stellar mass loss and mass transfer by taking into account the 
conservation laws of linear and angular momentum with respect 
to the centre of mass of the actual binary system. We do not 
make use of the orbital averages a priori, but give the variations 
as a function of the orbital phase, allowing different rates of mass 
loss and mass transfer along the orbit. In ^we calculate the rate 
of change of the eccentricity according to different assumptions 
of mass loss and accretion. In ^we evaluate the competition 
between the tidal circularisation and the eccentricity pumping 
due to our proposed mass loss, and show that the latter can prove 
effective in counteracting the former for systems such as Sirius 
and barium stars in which their AGB component did not fill its 
Roche lobe. A summary is presented in ^ 

2. Variational analysis of orbital parameters due to 
wind mass transfer or mass loss 

The dynamics of a binary orbit are well described by the two 
masses. Mi and M2, its energy, Eg^b, and its angular momentum, 
7oib (e g-, lGoldsteinill980l) . The fractional variation of the semi- 
major axis (also called the separation), a, and of the eccentricity, 
e, are expressed as 



6E, 



orb 



and 
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where we use the notation for the reduced mass problem, such 
that the total mass M = Mj + M2 and the reduced mass m = 
^'^^^ . For a detailed derivation of these expressions and those 
that follow in this section see Appendix lAl 



2. 1 . Angular momentum and energy 

In a system where r is the instantaneous distance between the 
stars and w is the instantaneous angular velocity of the compo- 
nents with respect to a co-moving inertial frame with its origin 
in the centre of mass (from here on the CM frame), the fractional 
variations of orbital angular momentum and energy are given by 
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6Eorb Sin 2a 6M 2a 6r 
Eoib m r M r r 



GM 



[2Vr6Vr + 2vg6vg] (4) 



respectively. We denote the velocity of one star relative to its 
companion by v = vi - V2, where v,- is the velocity of star / with 
respect to the CM. v also corresponds to the orbital velocity in 
the reduced mass problem. We also use the subscripts r and 6 to 
indicate the radial and transverse velocity components respec- 
tively. 

We consider a system in which both stars may lose mass in 
isotropic winds, part of which may be accreted by the companion 
such that during a time interval 6t an amount of mass 5Mi w is 
lost by star 1 and an amount 6M2,acc is accreted by star 2. In 
this way an amount of mass 5Mi_lost - ^A#i,w - 5Af2,ACC is 
effectively lost by star 2 from the system. To determine how the 
orbital velocity components vary we assume that star 2 receives 
an impulse when it accretes an amount of mass (5M2,acc coming 
from the wind of star 1, which at the moment it is accreted has a 
velocity wu in the CM frame. The corresponding expressions for 
mass lost from star 2 are obtained by interchanging the indices 1 
and 2. This yields the variations 
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where r,- coiTesponds to the distances of star / to the CM. The 
notation +(1 o 2) used here and below means: add terms to 
the left again with indices 1 and 2 interchanged. For quantities 
defined with an index in the left-hand side of the equation it is 
implicit that the corresponding expression for the other star is 
obtained by exchanging the indices in the right-hand side. By 
means of Eqs. (O, (|4]i, Q and (|6]l we can now write the orbital 
angular momentum and energy variations as 
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+ (1 o 2). 



In these equations vfrei.z = (wi2 - V2) is the relative velocity of 
the wind with respect to star 2 when accreted. We have also used 
the fact that there must be no displacements ^r,, which implies 
that 6rlr - (see AppendixlAlfor details). 

2.2. Behaviour of the wind 

To apply these expressions in calculations of orbital evolution, 
the detailed behaviour of the wind as it travels through the sys- 
tem must be followed from its point of origin until it is accreted 
or escapes the binary potential we ll. In general this re quires ei- 
ther ballistic calculations (e.g., iBrookshaw & Tavanilll993b or 



full-scale hydrody namical calculations (e.g.. lTheuns et al.lll996l : 
iNagae et al]|2004l) . For the purpose of this paper we assume the 
simple model described below. 

The velocity of the wind from star j when it is accreted by 
star / has a radial component, Wrji, that we will assume to be the 
velocity of the wind emanating isotropically from star /' as if it 
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were a single star, wj. For the transverse component, wgji, we 
assume that the specific angular momentum of the wind with re- 
spect to the CM is conserved and we express it as /iw,/ = Xj''/^- 
Note that the wind carries both the orbital and spin angular mo- 
mentum of the mass-losing star With the above definition = 1 
represents a point mass star, given that such an idealised object 
can only carry orbital angular momentum and no intrinsic spin. 
Once the wind is accreted by the companion star / a fraction 
of /iw,/ is transfered to the spin angular momentum of the accre- 

tor and the rest is transfered to the orbit. Thus wgjj - j^- ^ , 
where hor\, is the orbital specific angular momentum given by 
ho^^ - r'^ap- - GMa(l - e^) and = xi(l - Mi), which yields 



(wg,i2 - rioi) - rcj 
and similar for we.ai - r\(jj. 



M,2 



MiM 
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3. Rate of eccentricity change 

More useful than arbitrary variations of the physical parameters 
are their rates of change in time, thus from this point on we will 
transform 6 — > d/df. The distance and orbital velocity compo- 
nents are 
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where 6 is the orbital phase angle measured from the periastron. 
It is clear from its dependence on the orbital phase-angle that e 
is not constant along the orbit. This is true for constant mass loss 
and accretion rates, but especially if mass loss and accretion rates 
also depend on the orbital phase-angle. Employing the results we 
have derived so far we calculate using Eq. ^ 
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We note that the first term only depends on the assump t ion of 
isotropic mass loss and it is the same as given bv' Sokeil(l2000l) 
(see also Eggleton 2006). It shows that the eccentricity can in- 
crease if the mass-loss rate is higher near periastron than near 
apastron. The second term also depends on the additional as- 
sumptions about the wind behaviour (see ^2.21 ). These are com- 
plicated expressions that must be followed along the orbit if they 
are to have a useful meaning. The orbital period is given by 

P - -s^^U^TW) defining the orbital average of a pa- 

rameter A as (A) = J jf/'^df - J A^dO we compute the net 
rate of change of a and e in one orbital period. When an orbital 
average is applied to Eq. (fTTT i we obtain an expression which 



depends on (m,,lost), (m,;lost cos o), (m,- acc), (m,- acc cos 0), 
(M/,acc sin 0), (m,;acc cos^ 0) and (M/,acc sin^ 6*) for / = 1,2. 
These averages can be calculated only after M,lost(^) and 
Mi^AcciO) are known. 

3. 1 . Constant wind and constant mass accretion 

Usually the simplest case is assumed which corresponds to when 
MiACC = {Mi,ACc(0) ) and M;.lost = (M ;,LOST(g)) are constant 
along the orbit (e.g., iHurlev et alj|2002|) . This yields 



<£> _ A^2,ACC 

e M2 



'( I - Vl-e2)(l-e^)^iM 2 
e2 MiM 



2\ 



-I- (1 o 2). 



(12) 

Clearly with the assumptions about the wind made here the ec- 
centricity can only decrease, therefore if this effect is compara- 
ble to that of tidal friction it will only help to circularise t he orbit 
on a fa ster timescale, otherwise it can be neglected (.Hurlev et al.l 



3.2. Accretion from a fast isotropic wind 

We model the mass accretion with the iBondi & Hovia (11944 
mechanism, in which 

(GMifMi; 



M2,ACciO) = a- 



i,w 

'"A^Wrel.aVi 



(13) 



where Mi,w is the wind mass loss rate from star 1, the dis- 
tance to the source of the wind at the moment it was emit- 
ted and a is a numerical factor which parameterises the ef- 
ficiency of the accretion mechanism. In the case of a steady 
isotropic wind much faster than the orbital motion the trans- 
verse component of the wind velocity can be neglected, there- 
fore one can safely assume that Wiei,2 = wi and - r. Then 
the accretion rate Mi.ACcifi) - ■^2,oacc(1 + ecos0)^, where 

<M2.ACC(9)) 

((1+f cosfl)-) ■ 

the system is given by Mi lost - Mi,w - M2,acc- With this as- 
sumption the averaged rate of change of e is 



Mif) ACC = ^v^-'^-ALu^;/^ ^ rpj^g rnass from star 1 which is lost from 



<£) _ M2,() ACC ^ 

e M2 



2Mi(Mi +M2) 



[^jM2^ - 2Mi(M2 -H 2Mi)) -H (1 o 2). 



(14) 



In contrast to Eq.(fT2b. where e always decreases with time, 
Eq.(fT4l) shows that the eccentricity will increase if M2/M1 > 
(1 -I- VI + 12^1) /3fi (e.g.. Ml < 0.65M2 for the case of point 
masses, where ^1 = 1). 

3.3. Accretion from tfie slow wind of an AGB star 

The case of an AGB star with an accreting companion is some- 
what different from the two cases we have already reviewed. The 
velocity of the wind from an AGB star is of the order of the stel- 
lar escape velocity, which is also of the order of the orbital ve- 
locity, if the size of the AGB star is a significant fraction of its 



' iHurlev et alj ( l2002h have a different expression for e/e due to dif- 
ferent assumptions made in their calculations, however, their expression 
and ours are of the same order of magnitude. 
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Roche lobe. This implies that wi ~ -yj aa^e-) ^ ^° transverse 
velocity of the wind cannot be neglected as in the fast- wind case. 
Thus, making use of Eqs. (|9]l and ( fTOl i we write 



Wrel.2'(0) = 
Wl 



GM 



\2 



GM 



MiesinS 
M 



[^iM2^-Mi^](l -l-ecose) 
MiM 
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and by means of Eq. ( IT3T l we express the accretion rate as 

M2,()ACC(1 +e COS 0)2 
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Where M,oAcc = (^f (m,acc(^)) 
case. The orbital average can be only calculated numerically, but 



/ [fiM2^-Mi^](l+gcos9) y 
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(16) 
in this 



we find that when wi > 3 



CM 



the result of the integration is 

well approximated by the fast- wind assumption. The fact that the 
orbital phase can change significantly between the moment the 
wind was emitted and when it is accreted also introduces a phase 
angle shift in the relation between Ma.ACC and Mi,w. However, 
we performed numerical calculations of (e) /e by applying dif- 
ferent phase angle shifts in Eq. (fT3T l and the results do not differ 
by more than 15% from the zero phase-angle case. 



3.4. Mass loss from an AGB star 

A non-AGB star has a steep density gradient in its photospheric 
layers, i.e., it has a well defined radius. When the star is smaller 
than its Roche lobe any mass loss is in the form of a wind 
which may be enhanced by the presence of the companion (e.g., 
[Tout & Eggleton 1988). If such a star becomes larger than its 
Roche lobe then the mass loss rate is abruptly enhanced by or- 
ders of magnitude, governed by the Roche lobe overflow mecha- 
nism. Thus there is an almost discontinuous transition from one 
regime to the other and they can be treated as separate cases. 
AGB stars, on the other ha nd, have a shallow surf ace density 
gradient (e.g., Bowen 1988- |Pastetter & Ritterlll989) and there- 
fore the transition between the wind and the Roche lob e overflow 
mass-loss regimes is smooth. iKarovska et al.l (12005 ^ made ob- 
servations in X-rays of the symbiotic binary system Mira which 
show the AGB star is surrounded by material in the shape of 
a lobe as predicted by the Roche model. This suggests that the 
star is undergoing a transitional form of mass loss that can be 
described as wind Roche-lobe overflow. This provides support 
to the idea that the extended and weakly-bound atmosphere of 
an AGB star can b e highly influenced by the tida l force exerted 
by the companion. iFrankowski & Tvlendal ( 1200 Ih computed the 
enhancement of the wind mass-loss rate of giant stars in binary 
systems through the effect of the Roche potential on the local 
surface temperature and gravity. They find that the mass-loss 
enhancement depends on (R>f/Ri)^, where and Rl are ra- 
dius of the AGB star and its Roche-lobe radius, respectively. 
However, the overall effect is modest, less than a factor of 2, 
unless is very close to Their model does not account for 
the wind Roche-lobe overflow transition expec ted for AGB stars 
and suggested by the observations of Mira. ITout & EggletonI 



propose that the mass-loss rate of cool giants is enhanced 
by the presence of a companion, driven by the tidal friction 
torque. Thus, in their heuristic model the mass-loss enhance- 
ment depends on the sixth power of the Ri^ ratio. Building on 
this idea, we propose the following AGB star mass-loss model, 
which provides a smooth transition between the single-star wind 
mass-loss rate, Mj^ „, and the Roche-lobe overflow mass-loss 
rate, Mrlqf: 



R^ 

Rl 



(17) 



For eccentric orbits the Roche geometry does not apply so 
we utilise an instantaneous Roche-lobe which depends on the 
actual distance between the stars, r, as 



a(l 



(1 + ecosff) 



(18) 



where is the mass ratio of the AGB star to its companion and 

fiq) 

~ o.6?v3+in(i+^i/3)' following the approximation of|Eggleto^ 
(Il983h . With this prescription the mass loss rate along the or- 
bit is not constant, but enhanced when the system is at orbital 
phases close to the periastron. Thus, if the AGB star does not 
fill its Roche lobe in the periastron and there is no mass transfer 
(Mi^[) Acc - 0), the net contribution in one orbital period to the 
eccentricity variation rate is 



M \af{q^) 



R^ ^ 3(8 + \2e^ + e^) 
8(l-e2)''/2 ■ 



(19) 



Eq. ( fT9] l will always contribute to enhance the eccentricity, given 
that Mrlof is larger than Magb.w- The effect of mass transfer to 
the companion will in many cases act against the eccentricity en- 
hancement, but its contribution to Eq. ( fT9] l cannot be expressed 
analytically and it has to be calculated numerically (see ^3.2| and 
^3. 3b . However, it becomes significant only when a considerable 
amount of mass lost from the AGB star is transferred. 



4. Comparison of timescales 

In this section we calculate and compare the timescales on which 
the effects of circularisation and e ccentricity enhancem ent will 
take place, Tciic and respectivelv. lHurlev et al.l (l2002h give the 
following expression for the tidal evoluti on of e due to convec- 
tive damping by comb inin g the results of IZahnI d 19771 see also 
IZahn & Bouche3fT989l) andlHul ([1981 : 



18 /cc 



7 ^conv 



(\+q^) (R^ 
M» ^^2(1 _g2)13/2 

li(l-.2)3/2^4(.^)^ 

18 n 



(20) 



where Menv is the envelope mass of the AGB star and Tconv 
is the convective turnover timescale. A good approximation to 
the depth of the convective zone in an AGB star is its ra- 
dius, thus by employing Eq. (31) from iHurley et al.l (l2002h we 
obtain Tconv = 0.235 1 (Mgnv^^ ^/i^ )'^Vr, with U the lumi- 
nosity of the star and the parameters expressed in solar units, 
fljf s/n is the ratio of spin angular rotation rate of the AGB star, 
Q■l^ s, to the average orbital angular velocity, n. The functions 
fx(e^) are those derived by IHuti (119811) and the numerical fac- 
tor /conv = /'min(l,(Ptid/2Tconv)^) is from [Rasio et al.l d 19961) . 
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where Aid ' = |n - ^« M2n. IWrbunt & PhinnevI (Il995h found 
from an analysis of observed red-giant binaries that /' ^ 1, how- 
ever, this value may be different for the case of AGB stars. It is 
important to note that this calculation of Tciic gives a stronger cir- 
cularising effect for orbits with fini te eccentricit y values than the 
usually employed prescripti on by ZahnI ( 19771) which as sumes 
an almost circular orbit (e.g. jRasio et alJll996tlSokeiil200 0). 

We estimate the Roche-lobe o verflow mass-loss rate 
Mm.Q¥ following the presc r iptiori of RitteJ (r988l see als o 
iMever & Mever-Hofmeiste^ 119831; iPastetter & Ritte^ Il989h . 
where 



RLOF 



2n I fcfi 



■ill 



■p^,phF(^^-'), (21) 



with F(x) = 1 .23 -H 0.5 log X. The photospheric density of the 
AGB star is given by jKippenhahn & Weigertll990l;ISokej20"06l) 



P*,ph - 



2 niiifi^ GM^ 

3 ^Kf 



(22) 



where is the opacity of the envelope of the AGB star and 
and T>f are the mean molecular weight and effective temperature 
of the stellar surface, respectively, oth is the mass of the hydro- 
gen atom and and is the Boltzmann constant. A typical 1.5-2.5 
Mq AGB star has a surface temperature of about 3300K, a radius 
of about 300 R© and a luminosity of about 10000 L©. Most of the 
surface material in such a star is weakly ionised due to the low 
surface temperature, so we assume a typical surface molecular 
weight of 1.3 (see e.g., Soker 2006). According to the envelope 
opacity calculations of Marigoi (i2002). the corresponding opac- 
ity at this temperature is approximately - 2.5 x lO"'' cm^ g"' 
for densities of the order of 10"'^ - 10"" gem"-', which is 
the expected range of AGB surface densities. Employing this 
value of and the luminosity, radius, mass and mean molecu- 
lar weight mentioned above, the prescription of Eq. ^T\\ yields 
Mrlof ~ 10"^ Moyr"'. This is about 10000-30000 times larger 
than the single-star wind mass-loss rate of our typical star at the 
beginning of the AGB lifetime, calculated with the prescription 
of^assiliadis & Wood ( 1993). By the end of the AGB phase, at 
the onset of the superwind phase, Mrlof is still 50-200 times 
larger than the wind mass-loss rate. 

4. 1 . Eccentricity of the Sirius system 

Sirius A, a main s equence star of 2.1 Mp, and S irius B, a white 
dwarf of 1.05 Mq (iGatewood & Gatewoodlll978 !). form a binary 
system with period P = 50.1 yr and eccentricity e = 0.59 
jvan den BosI 11960 ^. Given the semi-major axis of the orbit, 
a = 20 AU, and the eccentricity of this system, the distance 
between its components at periastron is about 1700 R©. By the 
time the primary star (initially a star of about 5.5-6 M©) was in 
the AGB phase its radius was so large (R x 750 Rq) that it could 
have filled its Roche lobe during periastron passages. Whether 
or not the Roche lobe was filled, a strong tidal interaction must 
have occurred due to the proximity of the system components. 
We have tracked the evolution of binary systems which evolve 
into a main sequence star with the mass of Sirius A and a white 
dwarf with the mass of Sirius B, which from here on we will call 
Sirius-like systems. To compute the binary evolution models we 
employe d a version of the rapid synthetic code of binary evo- 
lution of Izzard et al. (2006), which makes use of modifications 
described in Bonacic Marinovic et al. (2007). The code employs 
Eq. |20] to model the effect of circularisation. The set of initial 
conditions for the models consists of the following: 




100 



a (AU) 



Fig. 1. Computed eccentricity of Sirius-like binary models with- 
out any eccentricity-enhancing mechanism as function of their 
separation, assuming the standard binary interaction. The filled 
circles indicate the orbital parameters of Sirius-like systems 
computed by assuming a tidal dissipation factor /' = 1. The 
open circles are similar, but assuming that the tidal dissipation is 
ten times weaker. The star indicates the observed orbital parame- 
ters of the Sirius system, which has a notably higher eccentricity 
than what standard models predict for its separation, even with 
substantially reduced tidal friction. 



- A primary with initial mass of 5.7 M©, which ends its life 
as a white dwarf similar to Sirius B, independently of the 
choice of mass loss. 

- 2 initial secondary masses (2.00 and 2.05 Mq) which yield a 
star similar to Sirius A. 

- 25 initial separations, oi, logarithmically separated between 
10 AU and 500 AU. 

- 25 initial eccentricities, separated linearly between 0.5 and 
0.99. 

- Solar metalHcity, Z = 0.02. 

In the standard scenario of tidal dissipation with constant 
mass loss along the orbit the eccentricity will always decrease. 
Fig.[T]shows the distribution of the synthetic Sirius-like systems 
in the a-e plane compared to the observed position of the Sirius 
system in this plane. The solid circles indicate the orbital prop- 
erties of Sirius-like systems when assuming f - 1, which show 
that the standard scenario of tidal interaction in binary evolu- 
tion predicts complete circularisation for systems separated up to 
about 30 AU. The closest system with an eccentricity of 0.59 has 
a separation of about 55 AU. The open circles show the orbital 
parameters predicted for Sirius-like systems when the strength 
of tidal dissipation is reduced to /' =0.1. The results are still 
far from reproducing the observed Sirius system. A separation 
of at least about 40 AU is required to obtain the observed eccen- 
tricity of 0.59. To reproduce the eccentricity and separation of 
Sirius within this scenario a tidal dissipation factor /' < 10"^ is 
needed, which is too strong a reduction to be acceptable. This 
demonstrates that the standard scenario of binary evolution can- 
not explain the observed eccentricity of the Sirius system, be- 
cause it predicts that the system should be circularised. 

However, this is not the case if the tidally enhanced mass- 
loss we propose in Eq. ( [TT] ) is considered. Due to its variable 
mass loss along the orbit an eccentricity enhancing effect coun- 
teracts the tidal circularisation. We compare Te and Tcirc to see 
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Fig. 2. Behaviour of eccentricity of a binary system with a 4.5M0 
AGB primary (initially a S.SM© star) and a 2.IM0 MS secondary 
as the system widens due to mass loss. The thick solid line in- 
dicates Te = Tciic when /' is assumed to be 0.5 and 5% of the 
mass lost via the enhanced mass-loss mechanism is accreted by 
the companion. The arrows indicate instantaneous directions in 
the diagram which a binary system will follow as it evolves. The 
star represents the current position in the a-e plane of the Sirius 
system. The thin dashed and dotted lines also indicate Te = Tdic, 
but considering /' = 0.25 and /' = 1, respectively. 
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Fig. 3. Computed eccentricity of Sirius-like binary models 
which take into account our proposed eccentricity enhancing 
mechanism plotted as function of their separation. The filled cir- 
cles indicate the orbital parameters of Sirius-like systems com- 
puted by assuming a tidal dissipation factor /' = 1 . The open 
circles and triangles indicate the results for similar calculations, 
but assuming that the tidal dissipation factor /' is 0.5 and 0.25, 
respectively. The star indicates the observed orbital parameters 
of the Sirius system, which are well reproduced if /' = 0.25 is 
assumed. 



which effect is dominant. We estimate their values at the mo- 
ment of evolution when the tidal dissipation effects are strongest, 
which is the bottleneck for ending up with an eccentric orbit. 
This occurs when the AGB star reaches its maximum radius, 700 
R0, at which point it has a luminosity of 38000 Lq and a core and 
envelope mass of 1.0 M© and 3.5 respectively. For the cal- 
culation of the Roche-lobe overflow mass-loss rate we assume 
ATy = 1 0~"^, which approximates the values calculated by lMarigol 
(I2OO2I) for the envelope opacity of high-metallicity stars with a 
temperature of about 3000 K. This gives Mrlof = 1-944 x 10"^ 
M0 yr'^ according to Eq. [21] while the wind mass-loss rate at 
this point is M» w - 5.168 x 10"^ Mq yr"'. While mass loss 
from a system tends to increase its separation, the tidal interac- 
tion tends to (pseudo-)synchronise the rotation of the AGB star 
with the orbit and thereby shrink the orbit due to transfer of or- 
bital angular momentum. This orbital shrinking effect is espe- 
cially strong in the beginning of the AGB phase, when the star 
expands rapidly and rotates slowly. However, once the AGB star 
rotation approaches a (pseudo-)synchronised rate with the orbit, 
the widening of the orbit due to mass loss dominates over the 
effect of tidal shrinking. Assuming this last regime, Fig.|2]shows 
how the a-e plane is divided into a region where the eccentricity 
decreases (upper-left area) and a region where the eccentricity is 
enhanced (lower-right area). Thus, if we adopt a modest reduc- 
tion of the tidal dissipation strength (/' - 0.5), any system with 
a separation larger than about 14 AU will eventually cross to the 
eccentricity-enhancing region, avoiding a circular orbit by the 
end of the AGB phase. In particular, the orbit of a system with 
orbital parameters similar to Sirius will remain eccentric. The 
thin dotted and dashed lines in Fig. |2] show that the location of 
the dividing line between decreasing and increasing eccentricity 
is sensitive to the tidal strength factor /' . This is because the ratio 
of timescales Tg/Tcirc only varies as {R^^laf-, i.e., a much weaker 
dependence than either of the individual timescales. There is 
a similar sensitivity to the adopted average opacity through 



^RLOF, which is also not certain, given that Eq. |2T| is only ap- 
proximate. Given these uncertainties, the results are not defini- 
tive. However, we have shown that by adopting reasonable val- 
ues for the uncertain parameters the orbit of Sirius can be kept 
significantly eccentric. Finally we note that the results are insen- 
sitive to the amount of mass accreted by the companion, as long 
as this is less than about 10-15% of the mass lost by the AGB 
star. 

So far we have shown how our eccentricity-enhancing mech- 
anism works by employing values of single-star models for the 
AGB primary at a given moment of evolution. However, the 
overall effect on the orbit of a given binary depends on how the 
relevant parameters change as a function of time as the binary 
evolves from zero age to its current state. We have therefore im- 
plemented the orbital evolution formul ae derived in Sect. [3 ] into 
the synthetic binary evolution code of llzzard etalJ (I2OOQ) . in- 
cluding our enhanced mass-l oss prescription. Full de tails of this 
implementation are given in jBonacic Marinovic & Pols 2QQ%. 
We calculate models for the same set of initial conditions as de- 
scribed above. The choice of the amount of mass that the com- 
panion accretes does not affect our results significantly, unless it 
is larger than about 15% of the total mass lost by the primary. 
We assume that 5% of the enhanced mass loss (Mj^ - M» ,„) is 
accreted by the companion, which increases the final mass of the 
secondary to that of Sirius A. Fig. [3] shows a comparison of the 
synthetic Sirius-like systems resulting from these models to the 
observed orbital properties of Sirius. The filled circles show the 
eccentricity and separation predicted by models which assume 
the standard tidal dissipation, f - I - Compared to the case when 
only tidal dissipation is considered (Fig. [U, the resulting bina- 
ries are significantly more eccentric in the separation range 15- 
40 AU, and the effect is comparable to reducing the tidal strength 
to /' = 0.1 in Fig.[T] The triangles in Fig. [3] show that to repro- 
duce the orbital properties of Sirius the tidal dissipation has to 
be weakened by only a factor of 4, which is still within a range 
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of reasonable values, as opposed to a reduction by more than a 
factor 1000 needed when only tidal circularisation is include. 

4.2. Eccentricities of fhe barium stars 

Barium stars show over-abundances of ^-process elements, most 
prominently of barium, and are found in binary systems with 
white dwarf companions. This constitutes evidence for mass ac- 
cretion from their companions when the latter were asymptotic 
giant branch (AGB) stars, given that i-process elements are syn- 
thesised in the AGB phase of evolution. Barium star binary sys- 
tems are observed to have periods between 80 and 10000 days 
and most of the systems with periods larger than about 600 days 
are significantly eccentric (see Fig. |4]i. These systems pose a 
problem to the standard binary evolution scenario which pre- 
dicts that orbits with periods shorter than about 3500-4000 days 
should have been circularised due to the tidal dissipation that 
must have taken place when the primary was an AGB star (e.g., 
iPols et af]|2003l) . The eccentricity enhancement resulting from 
our proposed AGB mass loss provides a mechanism which can 
effectively compete with tidal circularisation. We calculate the 
relevant timescales for a binary system with a 1.5 AGB pri- 
mary and a 1.1 secondary, with metallicity Z = 0.008. We 
assume that 10% of the mass lost by the AGB star is accreted 
by the companion, so that the latter accretes enough material to 
enhance its i-process element abundances significantly and be- 
come a barium star. The calculations are made at the AGB stage 
of evolution where the tidal dissipation is strongest (the circu- 
larisation bottleneck), where the radius of the AGB star is 280 
Rq, its luminosity is 7500 Lq and its core and envelope masses 
are 0.64 Mq and 0.6 respectively. Applying these values to 
Eq. |2T]we obtain that Mrlof - 1.75 x 10"-' Mq yr"' accord- 
ing to Eq. |2T1 while the wind mass-loss rate at this point is 
M^,w = 1.02X 10--**Mo3;r-'. 

Our results are shown in Fig . Bl and are compared to the 
observed data from iJorissen et alj (1 19981) in the P-e plane. The 
thick dashed line separates the region where the orbits will be 
circularised (to the left of the line) from the region where ec- 
centric orbits are possible (to the right of the line) in the stan- 
dard binary scenario, where only the tidal circularisation is act- 
ing. This limit is calculated by assuming r^jj^ = M/M and a 
tidal dissipation factor /' = 1, and reasonably well reproduces 
the li mit resul ting from binary population synthesis calculations 
(e.g., iPols e t al. 2003). The majority of the systems cannot be 
explained in this picture. The thick solid line indicates the same 
limit, but with the eccentricity pumping mechanism included, 
i.e., when ItJ ' -T^TiJ^,! - M/M. The evolution of binary systems in 
the P-e plane is similar to that in the a-e plane (shown in Fig.|2l), 
thus the eccentricities of the systems to the right of the thick line 
can be explained. This means that with the default tidal strength 
the eccentricity-enhancement mechanism cannot compete effec- 
tively with the tidal circularisation, except for a very limited pe- 
riod range around 3000 days. However, the tidal dissipation fac- 
tor /' is not very certain for AGB stars, given that it has been 
measured to be approximately unity m ostly using observations 
of red giants (IVerbunt & Phinnevlll99 5). Hence it is reasonable 
to calculate the circularisation limits for smaller values of /', as 
our calculations for the Sirius system suggest. Fig. |4] shows in 
thin lines the same limits that we have already described above, 
but assuming /' = 0.1. In this case the eccentricity of almost 
all systems which do not fill their Roche-lobes can be explained. 
[Karakas et al. (2000) found that with a similar choice of /' they 
could reproduce the eccentricities of all observed barium stars, 
including those with periods less than 1000 days. In their simu- 
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Fig. 4. Distribution in the P-e plane of barium star and the lim- 
its for circularisation set by different models. The cr osses indi- 
cate th e observational data of barium stars form Jorisse n et al.l 
(119981) . The thick lines delimit the regions where eccentric and 
circularised orbits are expected and are calculated for a system 
as described in the text, by assuming a tidal dissipation factor 
/' = 1 . The thin lines show the same limits as the thick lines, but 
this time calculated by assuming /' - 0.1. Solid lines indicate 
the circularisation limit when the eccentricity pumping mecha- 
nism we propose is considered. The dashed lines show the circu- 
larisation limit in the standard scenario. The dotted line indicates 
the limit for which Roche lobe overflow occurs at the periastron 
and the horizontal thin dotted Une indicates e - 0. 



lations they included a tidally enhanced AGB mass loss rate as 
well, but without any eccentricity-enhancing mechanism to com- 
pete against tidal circularisation. However, the results they report 
are difficult to understand in the light of our findings, especially 
because they excluded all systems in which the AGB star filled 
its Roche lobe. The systems to the left of the thick dotted line in 
Fig- SI fill their Roche lobes during the AGB and in most cases 
enter a common envelope (CE) phase. The evolution of the orbit 
during this phase is not a well understood problem and is also 
beyond the scope of this study. However, it seems more likely 
that an already eccentric system which enters a CE phase may 
remain eccentric than that a circular system becomes eccentric 
during the CE phase. Thus a mechanism which at least allows 
orbits to remain eccentric until the point when Roche-lobe over- 
flow occurs may be necessary to explain the eccentric systems 
which are to the left of the Roche-lobe limit. As shown in Fig. 
m assuming weak tidal friction with /' = 0. 1 makes it possible 
that systems with an AGB primary fill their Roche lobe with a 
significant eccentricity. 

It is important to note that the reduction of the tidal strength 
needed to explain the eccentricities of most observed barium- 
star systems which did not fill their Roche lobe (/' - 0.1) is 
calculated for one specific system, assumed to be representative. 
In order to test whether the eccentricity-enhancing mechanism 
explains the properties of the barium stars, the full evolutionary 
history has to be followed (as in the case of Sirius) for a vari- 
ety of system parameters in a binary population synthesis study. 
Furthermore, not only the orbital parameters have to be repro- 
duced, but also their abundance distributions. Such an analysis 
is beyond the scope of this work and we will present it in a subse- 
quent paper Preliminary results of population synthesis of bar- 
ium stars indicate that /' = 0.25 is needed for reproducing the 
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eccentricities of most observed systems which have not filled 
their Roche lobes, which is consistent with the reduction of tidal 
dissipation strength needed to reproduce the orbit of the Sirius 
system. 

5. Conclusions 

We have revised the description of the evolution of orbital pa- 
rameters due to mass loss and mass transfer in a binary system 
where the primary is an AGB star We also propose a tidally en- 
hanced mass-loss rate for AGB stars in binary systems which 
allows a smooth transition between wind mass loss and mass 
loss due to Roche lobe overflow. With our revised prescription 
for the evolution of the orbital eccentricity and the fact that our 
proposed mass loss is not constant along an eccentric orbit, we 
find an eccentricity enhancing mechanism which counteracts the 
circularising effect of tidal dissipation. We have shown that the 
standard scenario of binary interaction cannot explain the ob- 
served orbital parameters of the eccentric binary system Sirius 
unless the strength of tidal dissipation in AGB stars is at least 3 
orders of magnitude smaller than in normal red giants. On the 
other hand, our models with the eccentricity-enhancing mech- 
anism can reproduce the orbital properties of the Sirius sys- 
tem with a reasonable choice of the tidal dissipation parameter 
(/' - 0.25). Our eccentricity enhancing mechanism also allows 
binary systems containing barium stars to remain eccentric with 
periods as short as about 1000 days under reasonable parame- 
ter assumptions, while under the same assumptions in the stan- 
dard scenario of tidal circularisation only systems with periods 
longer than about 2500 days are expected to show a significant 
eccentricity. We also show that we can explain the eccentric- 
ities of almost all barium star systems which do not fill their 
Roche lobes if the uncertain convective tidal dissipation strength 
in AGB stars is reduced by a factor of 10, compared to what is 
measured for red giant stars. Moreover, this assumption allows 
Roche-lobe overflow to occur in significantly eccentric systems, 
which may be the key to explain shorter-period barium-star sys- 
tems which are still eccentric. 

Whether our eccentricity-enhancing mechanism is a vi- 
able solution to the problem of explaining short-period ec- 
centric binaries with (former) AGB stars, in particular barium 
stars, needs to be tested with a binary population synthesis 
model. Results of s uch a test will be given in a separate paper 
dBonacic Marinovrc & Pols.2007) . Preliminary barium star pop- 
ulation synthesis calculations indicate that when the full evo- 
lutionary history of a binary population is taken into account, 
a reduction of the tidal dissipation strength of only /' = 0.25 
is needed to reproduce the observed systems. This is consistent 
with our findings about the values of /' needed to reproduce the 
Sirius system. 

Acknowledgements. The authors would like to thank Frank Verbunt for pointing 
out that the problem of the eccentricity of the Sirius system is similar to that of 
the eccentricities of the barium stars. 
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Appendix A: Derivation of variation of orbital 
parameters 

The orbit of a binary star system can be described by the eccen- 
tricity e and the semi-major axis a. The change in these orbital 
parameters can be related to changes in the masses Mi 2 of the 
two stars, either due to mass loss or mass accretion, through the 
conservation laws for the total orbital energy Eo^b and the orbital 
angular momentum /orb, given by 



^orb 



GMm 
2a 



(A.l) 



Bonacic Marinovic et al.: Orbital eccentricities of binary systems with a former AGB star 



and 



(A.2) 



Here oj is the instantaneous angular velocity, r is the instanta- 
neous separation among the system components, M - M\ + M2 
is the total mass and m - M1M2/M is the reduced mass of the 
system. The variations of the orbital parameters a and e are then 
given by 

6a 6M 5m (J-Eorb 



M m Eoih 



and 



£orb ~ IT' ~' 



(A.3) 



(A.4) 



1 - 7orb -torb M 111 

In the following we will derive expressions for SEoib and dJorh- 
A.1. Variation of tlie Orbital Energy 

The total orbital energy is the sum of the kinetic energy and the 
gravitational energy. 



■^orb 



1 , 

—mv' 
2 



GM1M2 



(A.5) 



where r = |ri - ril is the instantaneous distance between the two 
stars and v = Vi - V2 is their relative velocity (see Fig. lA.lb . The 
variation of ( IA.5b is given by 



1 



SEorb - -mv — + 



Sm d(v^) 



GMm I5M 6m 6r 
r \ M m r 



). (A.6) 



The velocity v has components v, and vq - ruj in polar coor- 
dinates, so that the variation of 5(v^) can be written as 



^(v ) = 7.VrbVr H- Iru) (rSo) + a)6r) . 



(A.7) 



For later use, we will need to know the variation in terms of the 
variations for the individual components. 



6(v ) - Ivrdvr^i + 2r(jj(r\5(jJ + (jj5r\) H- (1 « 2). 



(A.8) 



The notation -i-(l o 2) means: add terms to the left again with 
indexes 1 and 2 interchanged. Note that this is an addition de- 
spite V being defined as the difference Vi -v^. The reason is that 
contrary to a Cartesian basis the direction of basis vectors in po- 
lar coordinates is not fixed but depends on the angular coordinate 

e. 

We will consider two effect s that change the mass of each 
star: Bondi-Hoyle accretion (Bo ndi & Hovlel 19441) and isotropic 
mass loss. The latter does not affect the velocity of the star rel- 
ative to the centre of mass, but the former does. Let 6M\ - 
SMi \cc denote the mass accreted by star 1 and let W21 be the 
velocity of the wind of star 2 at the moment it is accreted by star 
1 and W12 vice versa. Then conservation of momentum gives 

Midvi = (W21 - vi) 5Mi,Acc (A.9) 
and similar for star 2. Inserting this in ( |A.8t gives the relation 



2v ■ iW2l - Vi) H- (1 o 2). 

Ml 



(A. 10) 



It now remains to express 6M and 6m in terms of the variation 
of the individual masses 5M\ and 6M2- For the purpose of book- 
keeping we split the mass lost from star 1 (and similarly for star 



2) into two parts: the amount (5Mi_lost that is lost from the sys- 
tem and the amount SM2,acc that is eventually accreted by star 
2. Then 



(5M, 



-5Mi_L0ST - SM2,ACC + <JAf 1,ACC 



and 



SM2 = -(yM2,L0ST - 5Mi,ACC + 5M2,ACC- 



(A.ll) 



(A. 12) 



The variation of the total mass 5M is just 6M1 + 5M2 while the 
variation of the reduced mass is 



6m 
m 



5Mi,L0ST M2 



6M 



l.ACC 



l-^] + (l«2). 



Ml Ml + M2 Ml 

(A.13) 

If the transfer of momentum due to the wind is instantaneous the 
distance r between the two stars is not changed, hence 6r - 0. 
Inserting ( IA.8I ). ( lA.llI ). ( IA.12b and ( IA.13b into ( IA.6b then gives 



Eorb Ml 

5Mi_LOST 

Wi 

which is the same as ([8]). 
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(A. 14) 



A.2. Variation of the Orbital Angular Momentum 
The variation of the orbital angular momentum ( IA.2b is 
6Jorb 6m 6r 5vg 



^orb 



6m 6r 

h — 

m r 



(A. 15) 



where as before vg - rcj. Using ( IA.9I ), ( IA.13I ) and the fact that r 
does not depend on the variation of the two masses immediately 
gives 



5/orb SMit^cC 



J. 
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Ml 
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Ml we_2i 
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(A. 16) 
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which is the same as 



A.3. The Variation of the Semi-Major Axis 

The variation of the eccentricity that results from the above con- 
siderations is given in (fTTI) in the text. For completeness we 
here give the corresponding expression for the variation in a. 
Inserting (IA.14b and (IA.13b in ( IA.3b we obtain 
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